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Preface 


The  report  herein  presents  results  of  a  modeling  study  on  Bluestone  Reser¬ 
voir,  WV.  The  model  (CE-QUAL-W2)  was  used  to  determine  the  effects  of 
increased  pool  elevation  and  hydropower  retrofitting  cm  in-pool  and  release 
temperature  and  dissolved  oxygen.  This  report  was  prepared  in  die  Environ¬ 
mental  Laboratory  (EL),  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  Vicksburg,  MS.  The  study  was  sponsored  by  the  U.S.  Army  Engineer 
District,  Huntington,  and  was  funded  under  the  Military  Interdepartmental 
Purchase  Request  No.  E8593HW01  dated  27  October  1992. 

The  Principal  Investigators  of  this  study  were  Ms.  Dorothy  H.  Tillman  and 
Mr.  Thomas  M.  Cole  of  the  Water  Quality  and  Contaminant  Modeling  Branch 
(WQCMB),  Environmental  Processes  and  Effects  Division  (EPED),  EL.  This 
report  was  prepared  by  Ms.  Tillman  and  Mr.  Cole  under  the  direct  supervision 
of  Dr.  Mark  Dortch,  Chief,  WQCMB,  and  under  the  general  supervision  of 
Mr.  Donald  L.  Robey,  Chief,  EPED,  and  Dr.  John  Harrison,  Director,  EL. 
Technical  reviews  by  Drs.  Dortch  and  Bany  Bunch,  WQCMB,  are  gratefully 
acknowledged. 

At  the  time  of  publication  of  this  report,  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 

This  report  should  be  cited  as  follows: 

Tillman,  D.  H.,  and  Cole,  T.  M.  (1994).  “Bluestone  Phase  2 
temperature  and  dissolved  oxygen  modeling  study,”  Miscellaneous 
Paper  EL-94-2,  U.S.  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS. 


Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

acres 

4,046.873 

square  meters 

cubic  feet 

0.02831685 

cubic  meters 

feet 

0.3048 

meters 

mites  (U.S.  statute) 

1.609347 

kilometers 

square  mites 

9 eaftnno 

t .  wOwWO 

square  kilometers 

1  Introduction 


Background 

The  U.S.  Army  Engineer  District,  Huntington,  is  presently  considering 
raising  the  pool  at  Bluestone  reservoir  1 1  ft1  and  adding  conventional,  base¬ 
load  hydropower.  Through  the  Water  Operations  Technical  Support  (WOTS) 
program,  the  Huntington  District  contacted  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  Environmental  Laboratory  for  recommendations  on 
evaluating  effects  to  water  quality  if  the  proposed  modifications  were  made. 

Personnel  from  WES  met  with  the  Huntington  District  and  discussed 
evaluation  of  future  water  quality  conditions  at  Bluestone  Lake,  WV.  Recom¬ 
mendations  were  made  on  the  approach  to  determine  effects  of  project  modifi¬ 
cations  and  included  three  phases:  (a)  apply  the  SELECT  model  (Davis 
et  al.  1987)  to  evaluate  potential  dissolved  oxygen  (DO)  of  release  water  with 
hydropower  assuming  no  change  in  in-pool  conditions,  (b)  apply  the  time- 
varying,  two-dimensional  (laterally  averaged)  hydrodynamic  and  water  quality 
model,  CE-QUAL-W2,  to  evaluate  potential  changes  in  in-pool  and  release 
temperature  and  DO  assuming  a  gross  water  column  oxygen  demand  for  DO, 
and  (c)  apply  CE-QUAL-W2  with  all  water  quality  state  variables  activated  to 
more  accurately  define  potential  changes  in  future  in-pool  and  release  DO 
instead  of  having  to  make  broad  assumptions  about  the  depletion  rate. 

Phase  1  of  the  Bluestone  Water  Quality  Study  was  completed  by  personnel 
at  the  Huntington  District  with  guidance  from  the  WES  Hydraulics  Laboratory. 
The  WES  Environmental  Laboratory  conducted  Phase  2  as  requested  by  the 
Huntington  District  Results  from  Phase  2  are  presented  in  this  report 


Study  Objective 

The  Environmental  Laboratory  assisted  the  Huntington  District  by 
conducting  the  Phase  2  numerical  modeling  of  temperature  and  DO  in 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  vi. 
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Bluestone  Lake,  WV.  Model  results  from  CE-QUAL-W2  scenario  runs  were 
used  to  evaluate  potential  changes  in  in-pool  and  release  temperature  and  DO 
by  raising  the  pool  1 1  ft  and  adding  hydropower  to  the  project 


General  Modeling  Approach 

This  study  involved  applying  the  two-dimensional  (laterally  averaged) 
hydrodynamic  and  water  quality  model,  CE-QUAL-W2,  to  Bluestone  Lake  for 
temperature  and  DO  only.  DO  was  modeled  in  a  simplified  manner  using  a 
gross  water  column  oxygen  demand  (WOOD)  and  a  sediment  oxygen  demand 
(SOD).  This  approach  results  in  more  uncertainty  for  DO  predictions.  The 
DO  rate  parameters  were  adjusted  to  match  2  years  (a  wet  year.1983,  and  dry 
year,  1981)  of  observed  data.  The  assumption  in  this  approach  is  that  the 
change  in  pool  will  not  affect  the  WOOD  and  SOD  rates.  This  assumption  can 
not  be  confirmed  without  proceeding  to  the  recommended  third  phase.  The 
benefit  of  this  study  was  to  have  more  confidence  and  greater  resolution  (in 
terms  of  time  discretization  and  accuracy  of  release  DO  results)  than  the  first 
phase  study  recommended  by  WES  in  determining  impacts.  Sensitivity  analy¬ 
ses  were  also  run  by  adjusting  die  SOD  and  WOOD  rates  in  the  calibration 
and  verification  control  data  sets  to  see  which  parameter  had  a  greater  effect 
on  DO. 

After  calibration/verification,  two  scenario  runs  were  made:  (a)  raising  the 
pool  1 1  ft  and  (b)  raising  the  pool  1 1  ft  and  adding  hydropower.  Comparisons 
were  made  between  calibration/verification  results  and  scenario  results  for  both 
years  to  determine  impacts  to  temperature  and  DO  on  in-pool  and  release 
concentrations. 


Site  Description 

Bluestone  Dam  has  impounded  the  New  River  near  Hinton,  WV  (Figure  1), 
since  December  1949.  It  was  constructed  for  various  purposes,  including  flood 
control,  recreation,  and  fish  and  wildlife  enhancement  Two  major  tributaries 
drain  into  Bluestone  Lake,  New  River  and  Bluestone  River,  for  a  total  drainage 
area  of  4,565  square  miles.  At  normal  summer  pool  (1,410  ft  from  1  April 
through  29  November),  die  surface  area  of  2,039  acres  is  created  with  a  back¬ 
water  of  10.8  miles.  At  normal  winter  pool  (1,406  ft  from  1  December 
through  29  March),  the  surface  area  of  1,800  acres  is  created  with  a  backwater 
of  9.5  miles.  The  maximum  pool  elevation  for  flood  control  is  1,520  ft  and 
creates  a  backwater  of  36  miles.  A  mean  hydraulic  retention  time  of  6  days  is 
estimated  using  the  1985  growing  season  discharge  of  3,183  cfs. 

Bluestone  Dam  is  a  concrete  gravity  dam  structure  having  an  overall 
height  of  165  ft  with  the  top  elevation  at  1,535  ft  and  bottom  elevation  at 
1,369  ft.  Maximum  depth  of  the  reservoir  is  approximately  60  ft  for  normal 
summer  pool.  Discharge  is  through  16  gated  sluices  that  each  measure 
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Figure  1.  Biuestone  Lake 
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5  ft  8  in.  wide  by  10  ft  high,  and  the  center  line  of  the  openings  are  at  elevation  1,394  ft. 
Penstocks  have  been  installed  for  future  installation  of  hydropower  with  the  center  line  of 
intakes  at  elevation  1,383  ft  Discharge  from  Bluestone  Lake  ranges  from  approximately  1,000 
to  40,000  cfs. 
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Model  Description 


Model  Discussion 


CE-QUAL-W2  is  a  two-dimensional  model  that  predicts  vertical  and  longi¬ 
tudinal  variations  in  hydrodynamics,  temperature,  and  constituents  in  a  water 
body  through  time.  The  model  is  based  upon  the  Generalized  Longitudinal- 
Vertical  Hydrodynamics  and  Transport  (GLVHT)  model  of  rivers,  reservoirs, 
and  estuaries  (Buchak  and  Edinger  1984).  Earlier  versions  were  known  as  the 
Laterally  Averaged  Reservoir  Model  (LARM)  (Buchak  and  Edinger  1982). 
Development  of  the  GLVHT  model  has  been  ongoing  since  1975  by  WES  and 
J.  E.  Edinger  and  Associates  of  Wayne,  PA.  The  GLVHT  has  been  previously 
used  to  simulate  temperature  distributions  and  circulation  patterns  in  water 
bodies  and  has  been  applied  to  a  variety  of  systems  (Budiak  and  Edinger 
1984).  The  main  modifications  to  the  GLVHT  model  resulting  in  CE-QUAL- 
W2  were  the  inclusion  of  the  algorithms  to  simulate  water  quality  constituents. 

CE-QUAL-W2  is  based  upon  a  finite  difference  solution  of  the  laterally 
averaged  equations  of  fluid  motion  including  the  following:  (a)  the  free  water 
surface,  (b)  hydrostatic  pressure,  (c)  horizontal  momentum,  (d)  continuity, 

(e)  constituent  transport,  and  (f)  an  equation  of  state  relating  density  and  con¬ 
stituents  including  temperature  and  solids  concentrations  (dissolved  and 
suspended).  By  solving  for  the  water  surface  elevation  implicitly,  the  restric¬ 
tive  Courant  surface  gravity  wave  criterion  is  removed,  allowing  simulation  of 
reasonable  time  frames  for  field  applications,  such  as  entire  stratification 
cycles.  An  explicit  scheme  is  then  used  to  transport  heat  and  chemical/ 
biological  constituents.  The  model  has  the  capability  of  including  head  or 
flow  boundary  conditions,  branches,  multiple  withdrawals,  and  other  features 
that  allow  its  application  to  a  variety  of  situations. 

Basic  features  of  CE-QUAL-W2  are  summarized  below: 

a.  Two-dimensional  (laterally  averaged)  simulations  of  temperatures,  con¬ 
stituents,  and  flow  fields. 

b.  Hydrodynamic  computations  influenced  by  variable  water  density 
cruised  by  temperature  and  dissolved  and  suspended  solids. 
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c.  Simulation  of  the  interactions  of  numerous  biological/chemical  factors 
influencing  water  quality. 

d.  Allowance  for  multiple  inflow  loadings  and  withdrawals  from  tributar¬ 
ies,  point  and  nonpoint  sources,  precipitation,  branch  inflows,  and  out¬ 
flows  from  a  dam. 

e.  Allowance  for  multiple  branches. 

/.  Allowance  for  ice  cover  computations. 

g.  Allowance  for  variable  time  steps. 

h.  Allowance  for  flow  or  head  boundary  conditions,  making  it  applicable 
for  reservoir  or  estuarine  modeling. 

i.  Simulation  of  circulation  patterns. 

j.  Restart  capability. 

k.  Inclusion  of  evaporation  in  water  balance. 

/.  Heat  transfer  computations. 

m.  Variety  of  output  options. 

n.  Selective  withdrawal  capabilities. 

CE-QUAL-W2  conceptualizes  die  reservoir  as  a  grid  consisting  of  a  series 
of  vertical  columns  (segments)  and  horizontal  rows  (layers),  with  the  number 
of  cells  equal  to  the  number  of  segments  times  the  number  of  rows.  The  basic 
parameters  used  to  define  die  grid  are  the  longitudinal  spacing  (Ax,  in  meters) 
and  die  vertical  spacing  (h,  in  meters).  The  vertical  spacing  and  the  longitudi¬ 
nal  spacing  may  vary  spatially.  Each  cell  also  has  an  associated  width  that 
represents  an  average  value. 

CE-QUAL-W2  currently  simulates  20  water  quality  constituents  in  addition 
to  temperature  and  circulation  patterns.  Many  of  the  constituents  are  simulated 
simply  to  include  their  effects  upon  other  constituents  of  interest.  The  constit¬ 
uents  are  separated  into  four  levels  of  complexity,  permitting  flexibility  in 
model  application.  The  first  level  (Table  1)  includes  materials  that  are  conser¬ 
vative,  noninteractive,  or  do  not  affect  other  materials  in  the  first  level.  The 
second  level  (Table  1)  allows  the  user  to  simulate  the  interactive  dynamics  of 
oxygen-phytoplankton-nutrients.  The  third  level  (Table  1)  allows  simulation  of 
pH  and  carbonate  species,  and  the  fourth  level  allows  simulation  of  total  iron, 
which  is  important  during  anoxic  conditions.  The  model  calculates  in-pool 
water  volumes,  surface  elevations,  densities,  vertical  and  longitudinal  veloci¬ 
ties,  temperatures,  and  constituent  concentrations  as  well  as  downstream 
release  concentrations. 
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Table  1 

Water  Quality  Constituent  Levels 


1  Uvell 

Conservative  tracer 

Conform  bacteria 

Inorganic  suspended  solids 

Total  dissolved  solids  or  salinity 

Lav 

el  2 

Labile  dissolved  organic  matter 

Ammonia-nitrogen 

Refractory  dissolved  organic  matter 

Nitrate-nitrogen 

Phytoplankton 

Dissolved  oxygen 

Detritus 

Organic  sediments 

Phosphate-phosphorus 

Lav 

al  3 

Dissolved  inorganic  carbon 

Carbon  dnxide 

Alkalinity 

Bicarbonates 

pH 

Carbonates 

Laval  4 

Total  iron 

Data  Requirements 

CE-QUAL-W2  requires  a  database  that  includes  in-pool  initial  conditions, 
reservoir  geometry,  physical  coefficients,  biological  and  chemical  reaction 
rates,  and  time  sequences  of  hydrometeorological  and  inflowing  water  quality 
quantities.  Observed  release  water  quality  data  is  also  needed  to  evaluate  pre¬ 
dicted  release  conditions.  Calibration/verification  is  highly  dependent  on  the 
availability  of  in-pool  water  quality  constituent  concentrations  at  several  loca¬ 
tions  within  the  reservoir. 
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Calibration/Verification 


Calibration/Verification  Data  Sources 

The  model  was  calibrated  and  verified  for  a  dry  and  wet  water  year  (1981 
and  1983,  respectively).  The  different  data  types  necessary  to  calibrate  and 
verify  CE-QUAL-W2  for  the  Bluestone  system  were  as  follows: 

a.  In-pool  temperature  and  DO  data  for  various  stations  in  Bluestone 
Lake. 

b.  Release  data. 

c.  Bathymetry  data. 

d.  Tributary  inflow  rates  ami  constituent  concentrations. 

e.  Meteorological  data. 

/.  Water  surface  elevation  data. 

g.  Dam  outlet  specifications. 

h.  Reservoir  elevation-area-capacity  table. 

The  Huntington  District  provided  the  observed  in-pool,  release,  water  surface 
elevations,  and  calculated  inflow  data  for  die  2  study  years.  The  Huntington 
District  also  provided  the  sediment  range  survey  data  used  in  calculating  the 
reservoir  geometry,  elevation-volume  curve,  and  the  plans  from  the  proposed 
hydropower  study.  Inflow  temperature  data  were  obtained  from  CD  ROM  for 
the  U.S.  Geological  Survey  station  New  River  at  Glen  Lyn,  VA  (station 
number  03176500).  Meteorological  data  were  obtained  from  the  U.S.  Air 
Force  Environmental  Technical  Applications  Center  in  Asheville,  NC,  for  the 
Roanoke,  VA,  and  Beckley,  WV,  first-order  meteorological  stations.  Data 
from  the  Roanoke  station  were  used  for  calibration  and  verification  because  the 
Beckley  station  was  missing  data  for  the  year  selected  for  verification. 
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Observed  in-pool  data  were  available  on  a  monthly  basis  for  both  years. 
During  1981  (calibration),  observed  data  were  available  for  the  months  of 
April  through  September.  Consequently,  the  calibration  period  was  limited  to 
these  months.  Likewise  for  verification,  observed  data  were  available  only  for 
tire  months  of  May  through  October,  which  limited  the  simulation  period  to 
these  months. 


Calibration 

Before  actual  calibration  of  temperature  and  DO  could  be  conducted,  the 
water  balance  of  Bluestone  Lake  had  to  be  accomplished.  Adjustments  to  the 
bathymetry  data  and  the  elevation  of  the  bottom  datum  were  made  to  correct 
water  imbalances  in  the  system.  These  parameters  were  adjusted  until  the 
predicted  elevations  and  volumes  satisfactorily  matched  the  elevation-area- 
capacity  data  provided  by  the  Huntington  District  (Figure  2).  An  elevation- 
volume  relationship  was  also  developed  from  the  data  that  predicted  the  water 


ELEVATION  (M) 


Figure  2.  Elevation-volume  cun/e 

surface  elevations  (WSEL)  based  on  initial  reservoir  volume,  inflows,  and 
outflows  of  Bluestone  Lake.'  This  relationship  was  used  to  check  for  erroneous 
values  in  die  inflow  and  outflow  data.  Based  on  the  volume  change  resulting 
from  tire  values  of  inflows  and  outflows  being  used,  the  predicted  WSEL  for 
1981  did  not  match  what  had  been  observed.  The  predicted  WSEL  varied 
more  than  30  ft  for  some  days  where  the  measured  data  showed  very  little 
change.  After  consultation  with  personnel  at  the  Huntington  District,  it  was 
suggested  that  inflows  calculated  by  the  Huntington  District  be  used  instead  of 
the  U.S.  Geological  Survey  data  since  the  Glen  Lyn  station  had  problems 
during  the  1980s.  Once  the  calculated  inflow  data  were  used,  predicted  WSEL 
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were  well  within  the  0.5 -m  error  considered  acceptable  (Environmental  Labora¬ 
tory  and  Hydraulics  Laboratory  1986).  In  fact,  die  predicted  WSEL  were 
almost  a  perfect  overlay  of  the  observed  values  (Figure  3)  excluding  minor 
errors  (i.e.,  less  than  0.1  m  for  short  periods). 


JULIAN  DAY 


Figure  3.  Predicted  versus  observed  WSEL  for  1981 


Satisfactory  results  for  hydraulic  calibration  allowed  initiation  of  water 
quality  calibration.  Temperature  was  calibrated  first,  since  DO  is  temperature 
dependent  During  temperature  calibration,  adjustments  were  made  to  die 
Chezy  coefficient  and  wind  sheltering  coefficient  (Appendix  A).  They  were 
initially  set  to  values  recommended  in  the  user’s  manual  (Environmental  Labo¬ 
ratory  and  Hydraulics  Laboratory  1986).  Adjusting  these  parameters  improved 
temperature  predictions.  However,  only  after  restricting  the  lower  limit  of 
selective  withdrawal  to  elevation  1,387  ft,  was  the  thermocline  predicted  cor¬ 
rectly.  Bluestone  temperature  profiles  show  more  stratification  than  would  be 
expected  from  a  reservoir  having  such  a  short  retention  time  (approximately 
6  days  at  the  most).  For  instance,  hypolimnetic  temperatures  would  have  been 
expected  to  increase  as  the  summer  progressed.  However,  die  observed  data 
showed  very  small  changes  in  hypolimnetic  temperatures  throughout  the  sum¬ 
mer,  especially  in  1983  (see  Figures  4  and  5).  It  is  unclear  why  restriction  of 
selective  withdrawal  was  necessary,  but  it  was  originally  believed  that  a  coffer 
dam  was  in  place  upstream  of  the  dam.  After  checking  with  District  person¬ 
nel,  it  was  found  out  that  this  was  not  the  case.  Other  reasons  for  having  to 
restrict  the  selective  withdrawal  may  be  that  sedimentation  has  occurred  near 
the  dam  since  the  last  sediment  range  survey  or  groundwater  seepage  is 
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Figure  4.  Final  calibration  results  for  DO  and  temperature  (predicted  versus 
observed)  (Sheet  1  of  5) 
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cooling  the  hypolimnetic  waters.  District  personnel  may  want  to  investigate 
this  further. 

Temperature  and  DO  profile  results  for  calibration  are  shown  in  Figure  4. 
Ten  observed  temperature  and  DO  profile  data  stations  were  available  in  1981 
for  comparison  to  predicted  temperature  and  DO  profiles.  Table  2  lists  the 
observed  in-pool  stations  and  the  location  of  each  in  relation  to  Blues  tone 
Dam.  In  Figure  4,  DO  and  temperature  profiles  are  presented  for  each 
observed  Julian  day.  DO  profiles  are  plotted  first  beginning  with  the  most 
upstream  station  proceeding  in  the  downstream  direction  with  the  temperature 
profiles  for  the  same  day  plotted  next  in  the  same  order.  For  example,  in 
Figure  4,  the  first  observed  Julian  day  is  28  July  81,  and  die  order  of  the 
DO  profile  stations  is  1BLN20014,  1BLN20013,  1BLN20012,  1BLN20011, 
1BLN20010,  1BLN20009,  1 BLN 2000 3,  1BLN20008,  1BLN20007,  and 
1BLN20002.  The  temperature  profiles  for  that  day  follow  in  the  same  order. 


Table  2 

Observed  Profile  Stations 

Station  No. 

Distance 
from  Dam,  miles 

Segment  No. 

1981* 

1983* 

1BLN20002 

0.25 

31 

X 

X 

1 BLN20007 

1.00 

20 

X 

1BLN20008 

2.00 

28 

X 

X 

1BLN20OO3 

2.00 

26 

X 

X 

1 BLN20009 

4.00 

25 

X 

1BLN20010 

5.00 

22 

X 

X 

1 BLN2001 1 

6.00 

20 

X 

1BLN20012 

7.00 

18 

X 

X 

1BLN20013 

8.00 

17 

X 

1 BLN 200 14 

9.00 

15 

X 

X 

*  X  indicates  which  stations  were  available  for  that  year.  | 

Calibration  temperature  profile  predictions  for  all  stations  compared  favor¬ 
ably  with  tire  observed  data.  Initially,  inflow  temperature  boundary  conditions 
were  set  to  the  observed  Glen  Lyn  station  temperature  values.  Because  this 
station  was  approximately  15  miles  upstream  of  the  modeled  boundary  seg¬ 
ment,  tire  most  upstream  temperatures  were  being  overpredicted.  To  improve 
tire  upstream  temperature  predictions,  inflow  temperature  boundary  conditions 
were  set  to  observed  values  at  tire  most  upstream  station  (1BLN20014).  These 
values  were  more  realistic  to  use  as  boundary  conditions  and  helped  to 
improve  temperature  predictions  in  the  upper  reaches. 
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Although  DO  was  modeled  in  a  simplified  manner,  calibration  results  com¬ 
pared  favorably  with  observed  data  (Figure  4).  Since  there  were  no  observed 
inflow  DO  data  available  at  the  Glen  Lyn  station,  DO  boundary  conditions 
were  initially  assumed  to  be  saturated.  Using  saturated  DO  boundary  condi¬ 
tions  resulted  in  overprediction  of  the  most  upstream  DO.  DO  boundary  con¬ 
ditions  were  then  set  to  the  observed  values  at  the  most  upstream  station 
(1BLN20014).  Initial  DO  predictions  in  the  upper  reaches  were  improved, 
which  improved  DO  predictions  in  the  downstream  reaches  as  well. 

Further  calibration  of  DO  required  adjustments  to  the  SOD  and  WCOD 
rates.  Initially,  they  were  set  to  values  recommended  in  the  CE-QUAL-W2 
user’s  manual  (Environmental  Laboratory  and  Hydraulics  Laboratory  1986). 
The  SOD  and  WCOD  rates  were  not  varied  longitudinally,  but  were  set  the 
same  for  all  segments.  After  adjusting  the  SOD  and  WCOD  parameters,  DO 
profiles  were  improved  at  some  stations,  but  were  worse  at  others.  Since  there 
are  many  factors  (i.e„  inflow,  allochthonous  inputs,  algal  photosynthesis  and 
respiration,  and  wind)  influencing  DO  concentrations  throughout  a  reservoir 
(Cole  and  Hannan  1990),  it  was  decided  that  SOD  and  WCOD  rates  should  be 
varied  longitudinally.  DO  profile  predictions  were  then  significantly  improved 
throughout  the  reservoir.  Final  SOD  and  WCOD  rates  are  shown  in 
Appendix  A. 

Many  of  the  disparities  between  predicted  DO  and  observed  (especially  in 
the  epilimnion  on  25  July  1981  and  25  August  1981,  at  stations  1BLN20002, 
1BLN20007,  1BLN20008,  and  1BLN20003)  were  attributed  to  algal  produc¬ 
tion,  which  was  not  simulated  by  CE-QUAL-W2  during  this  phase  of  the 
study.  Since  DO  was  supersaturated,  the  higher  DO  values  observed  in  the 
epilimnion  could  not  be  predicted  without  the  inclusion  of  algae  as  a  modeled 
constituent 

Assessment  of  model  performance  for  release  conditions  was  conducted  by 
comparing  predicted  release  conditions  to  observed  conditions  at  a  station 
500  ft  downstream  of  the  dam.  Release  temperature  for  both  years  compared 
favorably  with  predicted  values;  however,  predicted  release  DO  was  consider¬ 
ably  lower  than  observed  values.  This  was  probably  due  in  part  to  reaeration 
and  the  inability  to  predict  the  higher  DO  values  in  tire  epilimnion  caused  by 
algal  production. 

Verification 

During  verification,  inflow  temperature  and  DO  boundary  conditions  were 
set  using  the  same  procedure  for  calibration.  All  other  parameters  (e.g.,  Chezy 
coefficient  and  wind  sheltering  coefficient)  were  also  set  the  same  as  during 
calibration  (Appendix  A).  This  included  having  to  restrict  selective  withdrawal 
at  the  same  elevation  to  correctly  predict  the  thermocline.  If  restricting  the 
selective  withdrawal  had  only  been  necessary  for  1  year,  that  doing  this  would 
have  been  suspect;  but  having  to  do  this  for  both  years  indicated  that  some¬ 
thing  is  influencing  the  temperature  profiles  in  this  region  of  the  reservoir. 
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Although  not  conclusive,  this  would  indicate  that  if  the  withdrawal  zone  was 
not  being  restricted,  anoxia  would  not  be  so  prevalent  in  this  zone.  Conse¬ 
quently,  the  system  would  be  flushed  out  with  such  a  short  residence  time. 

Results  for  verification  are  shown  in  Figure  S.  Six  observed  temperature 
and  DO  profile  data  stations  were  available  in  1983  for  comparison  of  predict¬ 
ed  values  (Figure  3).  Verification  results  are  plotted  the  same  as  calibration. 
Each  set  of  plots  has,  for  each  observed  Julian  day,  DO  profiles  plotted  first 
beginning  with  the  most  upstream  station  proceeding  in  the  downstream  direc¬ 
tion  with  the  temperature  profiles  plotted  next  in  the  same  order.  For  example, 
in  Figure  5,  for  the  first  observed  Julian  day,  the  order  of  the  DO  profile 
stations  is  1BLN20014,  1BLN20012,  1BLN20010,  1BLN20003,  1BLN20008, 
and  1BLN20002.  The  temperature  profiles  for  the  same  day  follow. 

An  acceptable  water  balance  was  obtained  for  verification.  The  predicted 
WSEL  was  well  within  0.5  m  tolerance  considered  acceptable  (Figure  6). 
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Figure  6.  Predicted  versus  observed  WSEL  for  1983 

Verification  temperature  and  DO  profile  predictions  for  all  stations  also  com¬ 
pared  favorably  with  the  observed  data  (Figure  S).  As  in  the  case  of  calibra¬ 
tion,  many  of  the  disparities  between  predicted  and  observed  DO  in  the 
epilimnion  (i.e.,  on  3  August  1983,  at  stations  1BLN20002,  1BLN20008, 
1BLN20003,  ami  1BLN20010)  were  attributed  to  algal  production.  Since 
sources  of  DO  other  than  reaeration  were  not  being  modeled,  the  higher  DO 
values  observed  in  the  epilimnion  could  not  be  predicted.  On  22  September 
1983,  DO  predictions  (Figure  5)  indicate  overturn  has  occurred  for  most  of  the 
reservoir  except  at  station  1BLN20002;  however,  this  was  not  indicated  by 
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observed  profile  data.  The  exact  date  of  overturn  is  difficult  to  predict  because 
of  limitations  in  meteorological  data  (i.e.,  met  stations  may  be  quite  a  distance 
from  the  project).  As  a  result,  predicted  overturn  is  often  a  few  days  off  from 
observed. 


Sensitivity  Analyses 

Sensitivity  analyses  were  performed  on  the  SOD  and  WOOD  rates  to  assess 
the  uncertainty  of  these  parameters  on  results  and  conclusions.  SOD  and 
WCOD  rates  were  increased  and  decreased  SO  percent  using  the  calibration/ 
verification  control  data  files.  Comparisons  were  made  between  calibration/ 
verification  results  (Figures  4  and  5)  and  results  from  the  sensitivity  analyses 
(Figures  B1-B8). 

Increasing  and  decreasing  the  SOD  rate  for  both  years  (Figures  B1  and  B2 
for  1981  and  Figures  BS  and  B6  for  1983)  showed  very  little  change  in  the 
predicted  DO  when  compared  with  calibration/verification  results.  This  was 
also  seen  in  the  release  DO  results  (Figures  B9  and  BIO).  The  SOD  results 
overlay  the  calibration  and  verification  results.  Since  the  SOD  only  affects  the 
DO  concentrations  at  the  sediment- water  interface,  these  results  are  reasonable. 

Adjustments  to  WCOD  rates  for  both  years  affected  DO  more  than  adjust¬ 
ments  to  the  SOD  rates  as  demonstrated  in  Figures  B3  and  B4  for  1981,  and 
Figures  B7  and  B8  for  1983.  When  WCOD  rates  were  increased  for  both 
years,  the  DO  values  in  die  entire  water  column  were  decreased  vertically  as 
well  as  longitudinally.  Increasing  WCOD  rates  caused  release  DO  values  to  be 
less  than  calibration/verification  release  (Figures  B9  and  BIO).  When  WCOD 
rates  were  decreased  (Figures  B4  and  B8)  for  both  years,  DO  values  in  the 
entire  water  column  were  increased  vertically  and  longitudinally.  This  caused 
the  release  DO  values  to  be  higher  in  comparison  to  calibration/verification 
release  results  (Figures  B9  and  BIO). 

Results  from  the  sensitivity  analyses  showed  that  DO  in  the  model  is  most 
sensitive  to  values  specified  for  WCOD. 
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4  Scenario  Results 


Changes  in  in-pool  and  release  conditions  were  assessed  by  comparing 
scenario  results  to  calibration  and  verification  results.  Two  proposed  modifica¬ 
tions  to  die  Bluestone  project  were  simulated.  Scenario  1  (Figure  Cl)  con¬ 
sisted  of  raising  the  pool  1 1  ft  Scenario  2  consisted  of  raising  the  pool  lift 
as  well  as  adding  hydropower  (Figure  C2).  For  all  runs,  no  reaeration  through 
the  sluice  gates  or  penstocks  was  assumed  to  occur.  In  Scenario  1,  die  dis¬ 
charge,  location,  arid  dimensions  of  the  intake  structure  were  assumed  to  be  tire 
same  as  in  calibration/verification  runs.  In  Scenario  2,  the  discharge  remained 
the  same  as  the  other  runs,  but  the  location  and  dimensions  of  the  intake  struc¬ 
ture  were  changed  to  conform  with  the  proposed  project  plans  for  Bluestone 
hydropower.  Selective  withdrawal  remained  restricted  during  the  scenario  runs 
since  calibration/verification  runs  indicated  this  was  necessary  to  simulate  the 
system.  If  selective  withdrawal  was  not  being  restricted.  Scenario  2  results 
may  be  slightly  different  since  the  penstock  location  is  deeper  in  tire  reservoir 
than  the  sluice  gates.  This  should  cause  DO  release  concentrations  to  be  lower 
than  calibration/verification  results. 

Comparisons  of  release  temperature  and  DO  between  calibration  (1981), 
Scenario  1,  and  Scenario  2  are  presented  in  Figure  C3.  Differences  between 
calibration  release  results  and  Scenario  1  and  Scenario  2  release  results  are 
shown  in  Figure  C4.  In  Figure  C4,  differences  were  calculated  as  calibration 
temperature  or  DO  minus  Scenario  1  temperature  or  DO  (represented  by  the 
dotted  line),  and  as  calibration  temperature  or  DO  minus  Scenario  2  tempera¬ 
ture  or  DO  (represented  by  the  dashed  line).  Similar  comparison  plots  and 
difference  plots  for  the  verification  year  (1983)  are  shown  in  Figures  C5 
and  C6,  respectively. 

Temperature  profile  results  from  Scenario  1  for  both  years  (Figure  Cl) 
demonstrate  that  raising  die  pool  11  ft  causes  the  thermocline  to  be  shifted 
deeper  in  die  reservoir.  This  causes  the  release  temperatures  for  both  years  to 
be,  on  the  average,  cooler  than  calibration/verification  results  until  around 
Julian  day  220  (Figures  C4  and  C6).  Although  the  thermocline  is  deeper  in 
the  reservoir  than  during  calibration  and  verification,  with  the  higher  pool,  it  is 
at  a  higher  elevation  in  relation  to  the  outlet  resulting  in  cooler  water  being 
withdrawn.  After  Julian  day  220,  temperature  releases  were,  on  the  average, 
warmer  (maximum  difference  1.1  °C)  than  calibration/verification  results.  This 


Chapter  4  Scenario  Results 


was  especially  true  for  the  dry  year  (1981).  Comparison  of  release 
temperatures  in  Figures  C3  and  C5  shows  that  adding  hydropower  had  very 
little  effect  on  release  temperature  results.  The  mean  release  temperatures  for 
the  calibration,  Scenario  1,  and  Scenario  2  simulations  using  1981  input  data 
were  24.92. 25.01.  and  25.03  4C.  respectively,  and  the  mean  release  tempera¬ 
tures  for  the  verification.  Scenario  1,  and  Scenario  2  simulations  using  1983 
input  data  were  25.02.  24.86,  and  24.84  °C,  respectively. 

DO  profiles  for  both  years  (Figure  Cl)  for  Scenario  1  show  that  because  of 
die  deeper  thermocline,  higher  DO  values  occur  deeper  in  the  reservoir.  This 
is  also  seen  in  Scenario  2  results  (Figure  C2).  Differences  in  DO  results 
shown  in  Figures  C4  and  C6  indicate  that,  on  the  average,  lower  DO  values 
(maximum  difference  approximately  5  mg/I  for  1981  and  3  mg/I  for  1983) 
were  released  for  both  scenarios  in  compaiison  with  catibration/verification 
releases.  The  mean  release  DO  concentrations  for  calibration.  Scenario  1,  and 
Scenario  2  simulations  for  1981  were  5.43, 4.82,  and  4.87  mg/I,  respectively, 
and  die  mean  release  DO  concentrations  for  verification.  Scenario  1,  and  Sce¬ 
nario  2  simulations  for  1983  were  5.88, 5.42,  and  5.42  mg/t,  respectively. 
Lower  DO  values  were  due  to  more  of  the  hypolimnetic  DO  being  available 
for  withdrawal.  The  greatest  DO  difference  between  the  two  scenarios  (Fig¬ 
ures  C4  and  C6)  occurs  between  Julian  day  255  and  day  265  for  both  years. 
This  difference  may  have  resulted  from  the  timing  of  overturn  and  the  differ¬ 
ence  in  the  withdrawal  zone  caused  by  the  different  intake  locations  and 
dimensions. 
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5  Summary  and  Conclusions 


CE-QUAL-W2  was  applied  to  Bluestone  Lake,  WV,  to  evaluate  impacts  to 
in-pool  and  release  temperature  and  DO.  The  model  was  calibrated  and  veri¬ 
fied  for  a  dry  and  wet  year  (1981  and  1983,  respectively).  After  calibration/ 
verification,  sensitivity  analyses  were  perfoimed  on  tire  SOD  and  WOOD  rates. 
Two  scenario  runs  were  simulated  looking  at  (a)  raising  the  pool  1 1  ft  and 
(b)  raising  the  pool  1 1  ft  and  adding  hydropower. 

Raising  the  pool  11  it  and  adding  hydropower  caused  changes  in  both 
in-pool  and  release  temperature  and  DO  when  compared  with  calibration  and 
verification  results.  From  the  two  scenarios  simulated,  the  following  conclu¬ 
sions  were  derived: 

a.  Temperature  [nofiles  for  most  stations  (especially  stations  closer  to  the 
dam)  showed  deeper  thermoclines  resulting  in  higher  DO  values  deeper 
in  the  reservoir.  Release  temperatures  increased  as  much  as  1.1  °C. 
Most  of  the  higher  release  temperatures  occurred  during  die  latter  half 
of  the  simulation  for  both  years.  Average  release  temperatures  for  the 
simulation  period  were  similar  in  value  between  calibration/verification 
and  the  scenario  results. 

b.  The  average  decrease  in  DO  releases  was  approximately  0.6  mg/f  for 
both  years.  Decreases  in  release  DO  occurred  throughout  the  simula¬ 
tion  period. 

c.  The  addition  of  hydropower  (Scenario  2  Figure  C2)  did  not  signifi¬ 
cantly  affect  temperature  and  DO  results  when  compared  with  Sce¬ 
nario  1  results  (Figure  Cl).  Selective  withdrawal  was  restricted  for 
these  runs  since  this  was  necessary  to  calibrate  and  verify  the  model. 
Late  in  the  study,  Scenario  2  was  rerun  for  both  years  with  selective 
withdrawal  not  being  restricted.  In  these  runs,  temperature  profile 
results  for  both  years  showed  very  little  thermal  stratification.  In  addi¬ 
tion,  IX)  profile  concentrations  were  also  higher  and  deeper  in  the 
reservoir,  and  the  mean  release  DO  concentrations  were  about  the  same 
as  Scenario  2  results.  The  plots  from  these  runs  have  not  beat  included 
in  the  report,  but  can  be  obtained  upon  request. 
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Blueatone  Reservoir  Control  File  for  CE-QUAL-W2 

TITLE  C . TITLE . 

Bluestone  Reservoir  calibration  -  run  24 
Adjusted  inflow  temperatures  &  DO  concentrations 
Hyd  -  CHEZY  -  50.0,  WSC  -  0.7,  KBSW  -  28,  dltf  -  0.5 
Hq  -  sod  ~  variable,  WOOD  -  variable 


TIME  CON 

TMSTRT 

TMEND 

YEAR 

180.85 

265.85 

1981 

DLT  CON 

NOT 

MINDLT 

1 

1.0 

DLT  DATE 

DLTD 

DLTD 

DLTD 

DLTD 

0.0 

DLT  MAX 

DLTMAX 

3600.0 

DLTMAX 

DLTMAX 

DLTMAX 

DLT  FRN 

DLTF 

0.50 

DLTF 

DLTF 

DLTF 

SURFACE 

KT 

DATUM 

15 

417.3 

BRANCH  G 

US 

DS 

UHS 

DHS 

Br  1 

2 

31 

0 

0 

Br  2 

34 

37 

0 

29 

LOCATION 

LAT 

LONG 

37.6 

80.9 

INIT  CND 

IT2 

IICETH 

WTYPE 

-1.0 

0.0 

FRESH 

CALCULAT 

VBC 

MBC 

PQC 

PQTC 

OFF 

OFF 

OFF 

OFF 

INTERPOL 

INFIC 

TRIC 

DTRIC 

HDIC 

ON 

OFF 

OFF 

OFF 

DEAD  SEA 

WINDC 

QINC 

QOUTC 

HEATC 

ON 

ON 

ON 

ON 

ICE  COVER 

ICEC 

SLICE 

SIflTEX 

ALBEDO 

OFF 

DETAIL 

TERM 

0.25 

TRANSPORT 

SLTRC 

THETA 

QUICKEST 

0.00 

WSC  NUMB 

NWSC 

1 

WSC  DATE 

WSCD 

0.000 

WSCD 

WSCD 

WSCD 

WSC  COEF 

WSC 

0.70 

WSC 

WSC 

WSC 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

PHIO 

3.142 

3.142 

EVC 

OFF 

PRC 

OFF 

OUTIC 

ON 

WDIC 

OFF 

METIC 

ON 

HWI 

10.0 

BETAI 

0.6 

GAMMAI 

0.07 

ICEMIN 

0.05 

ICET2 

4.0 

WSCD 

WSCD 

WSCD 

WSCD 

WSCD 

WSC 

WSC 

WSC 

WSC 

WSC 

HYD  COEF 

AX 

1.0 

IOX 

1.0 

AZMIN 

1.4E-6 

DZNIN 

1.4E-7 

DZMAX 

1000.0 

CHEZY 

50.0 

SSL  WITH 

swc 

ON 

swc 

OFF 

SWC 

SWC 

SWC 

SWC 

SWC 

SWC 

SWC 

N  STRUC 

NSTR 

1 

NTSR 

1 

NSTR 

NSTR 

NSTR 

NSTR 

NSTR 

NSTR 

NSTR 

K  BOTTOM 

XBSW 

28 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

SINK  TYPE 
Br  1 

Br  2 

SINX 

LINE 

SINK 

SINK 

SINK 

SINK 

SINK 

SINK 

SINK 

SINK 

E  STRUC 

Br  1 

Br  2 

ESTR 

426.5 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

W  STRUC 

Br  1 

Br  2 

WSTR 

240.85 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

M  OUTLET 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

0  LAYER 

Br  1 

Br  2 

ROUT 

ROOT 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

N  WDRHAL 

NWD 

0 

W  SEGMNT 

IWD 

0 

IWD 

IWD 

IWD 

IWD 

IWD 

IWD 

IWD 

IWD 

W  LAYER 

KWD 

0 

KWD 

KWD 

KWD 

KWD 

KWD 

KWD 

KWD 

KWD 

H  TRIBS 

NTR 

0 

TRIB  SEG 

ITR 

29 

ITR 

ITR 

ITR 

ITR 

ITR 

ITR 

ITR 

ITR 

DST  TRIB 

DTRC 

OFF 

DTRC 

OFF 

DTRC 

DTRC 

DTRC 

DTRC 

DTRC 

DTRC 

DTRC 

SNAPSHOT 

FORM 

LONG 

UPRNC 

OFF 

WPRNC 

OFF 

TPRNC 

ON 

SHRT  SEG 

IPRSF 

2 

36 

IPRSF 

5 

37 

IPRSF 

10 

IPRSF 

15 

IPRSF 

20 

IPRSF 

25 

IPRSF 

30 

IPRSF 

31 

IPRSF 

35 

LONG  SEG 

XPELF 

2 

20 

XPRLF 

4 

22 

IPRLF 

6 

25 

IPRLF 

8 

26 

IPRLF 

10 

28 

IPRLF 

12 

29 

IPRLF 

15 

30 

IPRLF 

17 

31 

IPRLF 

18 
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SNP  PRUT 

SNPC 

ON 

NSNP 

4 

SNP  OATS 

SNPD 

180.85 

SNPD 

208.85 

SNPD 

236.85 

SNPD 

265.85 

SNP  FREQ 

SNPF 

100.0 

SNPF 

2.0 

SNPF 

2.0 

SNPF 

100.0 

FRF  PLOT 

PRFC 

ON 

NPRF 

4 

NIPRF 

10 

PUT  DATE 

PRFD 

180.85 

PRFD 

208.85 

PRFD 

236.85 

PRFD 

265.85 

PRF  FREQ 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

100.0 

PRF  SEG 

IPRF 

15 

31 

ZPRF 

17 

IPRF 

18 

ZPRF 

20 

TSR  PLOT 

TSRC 

ON 

NTSR 

1 

TSR  DATE 

TSRD 

181.5 

TSRD 

TSRD 

TSRD 

TSR  FREQ 

TSRF 

0.25 

TSRF 

TSRF 

TSRF 

VPL  PLOT 

VPLC 

OFF 

NVPL 

4 

VPL  DATE 

VPLD 

181.85 

VPLD 

209.85 

VPLD 

237.85 

VPLD 

266.85 

VPL  FREQ 

VPLF 

100.0 

VPU 

100.0 

VPU 

100.0 

VPU 

100.0 

CPL  PLOT 

CPLC 

OFF 

NCPL 

1 

CPL  DATE 

CPLD 

212.208 

CPLD 

CPLD 

CPLD 

CPL  FREQ 

CPU 

1.0 

CPU 

CPU 

CPU 

RESTART 

RSOC 

OFF 

NRSO 

1 

RSZC 

OFF 

RSO  DATE 

RSOD 

267.85 

RSOD 

RSOD 

RSOD 

RSO  FREQ 

RSOF 

100.0 

RSOF 

RSOF 

RSOF 

CST  COMP 

CCOHPC 

ON 

LZMC 

OFF 

SDC 

ON 

FREQUK 

1 

SNPD 

SNPD 

SNPD 

SNPD 

SNPD 

SNPF 

SNPF 

SNPF 

SNPF 

SNPF 

PRFD 

PRFD 

PRFD 

PRFD 

PRFD 

PRFF 

PRFF 

PRFF 

PRFF 

PRFF 

ZPRF 

22 

ZPRF 

25 

ZPRF 

26 

ZPRF 

28 

IPRF 

29 

TSRD 

TSRD 

TSRD 

TSRD 

TSRD 

TSRF 

TSRF 

TSRF 

TSRF 

TSRF 

VPLD 

VPLD 

VPLD 

VPLD 

VPLD 

VPU 

VPU 

VPU 

VPU 

VPU 

CPLD 

CPLD 

CPLD 

CPLD 

CPLD 

CPU 

CPU 

CPU 

CPU 

CPLF 

RSOD 

RSOD 

RSOD 

RSOD 

RSOD 

RSOF 

RSOF 

RSOF 

RSOF 

RSOF 
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CST  ACT 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CST  ICON 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-2.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CST  PRNT 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CIN  CON 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CTR  CON 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CDT  CON 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CFR  CON 

FRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

EX  COEF 

EXH20 

EXINOR 

EXORG 

BETA 

0.45 

0.01 

0.3 

0.45 

CO LI FORM 

COLQIO 

COLDK 

1.04 

1.4 

S  SOLIDS 

SSETL 

2.0 

ALGAE 

AGROW 

AMORT 

AEXCR 

ARESP 

ASETL 

ASATUR 

ALGDET 

1.5 

0.05 

0.02 

0.02 

0.14 

50.0 

0.80 

ALG  RATE 

AGT1 

AGT2 

AGT3 

AGT4 

AGK1 

AGK2 

AGK3 

AGK4 

10.0 

30.0 

35.0 

40.0 

0.1 

0.98 

0.98 

0.1 

DISS  ORG 

LABDK 

LRFDK 

REFDK 

0.12 

0.001 

0.001 

DETRITUS 

DETDK 

DSETL 

0.08 

0.35 

ORG  RATE 

OMT1 

OMT2 

OMK1 

0MK2 

4.0 

20.0 

0.1 

0.98 
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SEDIMENT  SEDDK 

0.10 


S  DEMAND 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

0.050 

0.050 

0.050 

0.10 

0.150 

0.20 

0.250 

0.250 

0.250 

0.250 

0.30 

0.350 

0.350 

0.350 

0.350 

0.350 

0.350 

0.30 

0.250 

0.250 

0.250 

0.250 

0.250 

0.20 

0.20 

0.20 

0.20 

0.150 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.150 

0.150 

0.10 

0.050 

WCOD  TEMP 

TWCOD 

1.0147 

HC  DEMAND 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.40 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

CBOD 

KBOD 

TBOD 

RBOD 

0.25 

1.0147 

1.85 

PHOSPHOR 

P04REL 

PARTP 

AHSP 

0.015 

1.2 

0.009 

AMMONIA 

NH3REL 

NH3DK 

PARTN 

AHSN 

0.08 

0.12 

1.0 

0.014 

NH3  RATE 

NH3T1 

NH3T2 

NH3K1 

NH3K2 

5.0 

20.0 

0.1 

0.98 

NITRATE 

N03DK 

0.12 

N03  RATE 

N03T1 

N03T2 

N03K1 

N03K2 

5.0 

20.0 

0.1 

0.98 

SED  C02 

C02REL 

0.1 

IRON 

FEREL 

FESETL 

0.5 

2.0 

STOICHMT 

02NH3 

020RG 

02RESP 

02ALG 

BIOP 

BION 

BIOC 

4.57 

1.4 

1.1 

1.4 

0.011 

0.08 

0.45 

02  UNIT  02LIM 

0.00 


BTH 

FILE . 

. BTHFN 

bth.npt 

VPR 

FILE . 

. VPRPH 

vpr.npt 

LPR 

FILE . 

. LPRPK 

lpr.npt 

RSI 

FILE . 
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rai.npt 

KPT1  . 

. METFN. 

aet.npt 

QWD  FILE . 

. . QWDFH  -  . 

not  used 

QIN  FILE . 

Br  1  qin_brl.npt 

Br  1  qin~br2.npt 

TIM  FILE . 

. TINFN. . 

Br  l  ti.n_brl.npt 

Br  l  tin_br2.npt 

CIN  FILE . 

. CINFN. . 

Br  1  cin_brl.npt 

Br  l  cin_br2.npt 

QOT  FILE . 

Br  1  qot_brl.npt 

Br  1  qot_br2.npt 

QTR  FILE . 

Tr  l  qtr_trl.npt 

TTR  FILE . 

Tr  i  ttr_trl.npt 

CTR  FILE . 

Tr  l  ctr_trl.npt 

QDT  FILE . 

. QDTFN. - 

Br  1  not  used 

Br  1  not  used 

TOT  FILE . 

Br  l  not  used 

Br  1  not  used 

COT  FILE . 

. CDTFN. . 

Br  l  not  used 

Br  1  not  used 

PRE  FILE . 

. PREFN. . 

Br  l  not  used 

Br  1  not  used 

TPR  FILE . 

. TPRFN- . 

Br  1  not  used 

Br  1  not  used 

CPR  FILE . 

. CPRFH. . 

Br  l  not  used 

Br  1  not  used 

EUH  PILE . 

. EtTOFN.  - 

Br  1  not  used 

Br  l  not  used 
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TUH  FILE . 

Br  i  not  uaad 

Br  1  not  used 

CUH  FILE . 

Br  1  not  ussd 

Br  1  not  ussd 

BOH  FILE . 

Br  1  not  ussd 

Br  1  not  ussd 

TOH  FILE . 

Br  1  not  ussd 

Br  1  not  ussd 

CDH  FILE . 

Br  I  not  ussd 

Br  1  not  ussd 

BNP  FILE . 

snp_run28.opt 

TSR  FILE . 

tsr_run28.opt 

PRF  FILE . 

prf_run28 . opt 

VPL  FILE . 

vpl.opt 

CPL  FILE . 

cpl.opt 


TUHFN 

CUHFN 

EDHFN 

TDHFN 

CDHFN 

SNPFN 

TSRFN 

PRFFN 

VPLFN 

CPLFN 
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Bluestone  Reservoir  Control  File  for  CE-QUAL-W2 

TITLE  C . TITLE . 

Blueetone  Reservoir  1983  verification  -  run  25 

Adjusted  water  surface,  inflow  temperatures,  t  DO  concentrations 
Hyd  -  CHEZY  -  50.0,  WSC  -  0.7,  KBSW  -  28,  dltf  -  0.9 
HQ  -  SOD  “  variable,  WCOD  -  variable 


TINE  CON 

TNSTRT 

TKEND 

YEAR 

158.8 

264.85 

1983 

DLT  CON 

NUT 

MXNDLT 

1 

1.0 

DLT  DATE 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

0.0 

DLT  MAX 

DLTMAX 

3600.0 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLT  FRN 

DLTF 

0.90 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

SURFACE 

KT 

DATUM 

15 

417.30 

BRANCH  6 

US 

DS 

UHS 

DHS 

PHIO 

Br  1 

2 

31 

0 

0 

3.142 

Br  2 

34 

37 

0 

29 

3.142 

LOCATION 

LAT 

LONG 

37.6 

80.9 

INIT  CND 

IT2 

IICETH 

HTYPE 

-1.0 

0.0 

FRESH 

CALCULAT 

VBC 

MBC 

PQC 

PQTC 

EVC 

PRC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

INTERPOL 

INFIC 

TRIC 

DTRIC 

HDIC 

OUTIC 

HDIC 

METIC 

ON 

OFF 

OFF 

OFF 

ON 

OFF 

ON 

DEAD  SEA 

WINDC 

QINC 

QOUTC 

HEATC 

ON 

ON 

ON 

ON 

ICE  COVER 

ICEC 

SLICE 

SLHTEX 

ALBEDO 

HHI 

BETAI 

GAMMAI 

XCEMIN 

ICET2 

OFF 

DETAIL 

TERM 

0.25 

10.0 

0.6 

0.07 

0.05 

4.0 

TRANSPORT 

SLTRC 

THETA 

QUICKEST 

0.00 

WSC  NUMB 

NWSC 

1 

HSC  DATE 

WSCD 

0.000 

HSCD 

HSCD 

HSCD 

HSCD 

HSCD 

HSCD 

HSCD 

WSCD 

HSC  COEF 

HSC 

0.70 

HSC 

HSC 

HSC 

HSC 

HSC 

HSC 

HSC 

WSC 
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KYD  COEF 

AX 

1.0 

IDX 

1.0 

AZMIN 

1.4E-6 

DZMIN 

1.4E-7 

DZMAX 

1000.0 

CHEZY 

50.0 

SEL  WITH 

swc 

ON 

SWC 

OFF 

SWC 

SWC 

SWC 

SWC 

SWC 

SWC 

SWC 

N  STRUC 

NSTR 

1 

NTSR 

1 

NSTR 

NSTR 

NSTR 

NSTR 

NSTR 

NSTR 

NSTR 

K  BOTTOM 

KB SW 
28 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

KBSW 

SINK  TYPE 
Br  1 

Br  2 

SINK 

LINE 

SINK 

SINK 

SINK 

SINK 

SINK 

SINK 

SINK 

SINK 

E  STRUC 

Br  1 

Br  2 

ESTR 

426.5 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

ESTR 

W  STRUC 

Br  1 

Br  2 

WSTR 

240.85 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

WSTR 

N  OUTLET 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

NOUT 

O  LAYER 

Br  1 

Br  2 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

KOUT 

N  WDRWAL 

NHD 

0 

W  SEGMNT 

IWD 

0 

IWD 

IWD 

IWD 

IWD 

IWD 

IWD 

IWD 

IWD 

W  LAYER 

KWD 

0 

KWD 

KWD 

KWD 

KWD 

KWD 

KWD 

KWD 

KWD 

N  TRIBS 

NTR 

0 

TRZB  SEG 

ITR 

29 

ITR 

ITR 

ITR 

ITR 

ITR 

ITR 

ITR 

ITR 

DST  TRIB 

DTRC 

OFF 

DTRC 

OFF 

DTRC 

DTRC 

DTRC 

DTRC 

DTRC 

DTRC 

DTRC 

SNAPSHOT 

FORM 

LONG 

UFRNC 

OFF 

WPRNC 

OFF 

TPRNC 

ON 

SHRT  SEG 

IPRSF 

2 

36 

IPRSF 

5 

37 

IPRSF 

10 

IPRSF 

15 

IPRSF 

20 

IPRSF 

25 

IPRSF 

30 

IPRSF 

31 

IPRSF 

35 

LONG  SEG 

IPRLF 

2 

20 

IPRLF 

4 

22 

IPRLF 

6 

25 

IPRLF 

8 

26 

IPRLF 

10 

28 

IPRLF 

12 

29 

IPRLF 

15 

30 

IPRLF 

17 

31 

IPRLF 

18 
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SNP  PUNT 


SNPC 

ON 


NSNP 

7 


SNP  DATE 

SNPD 

144.80 

SNPD 

158.80 

SNPD 

172.80 

SNPD 

200.80 

SNPD 

214.80 

SNPD 

228.80 

SNPD 

264.80 

SNPD 

284.80 

SNPD 

SNP  FREQ 

SNPF 

100.0 

SNPF 

100.0 

SNPF 

100.0 

SNPF 

100.0 

SNPF 

100.0 

SNPF 

100.0 

SNPF 

100.0 

SNPF 

100.0 

SNPF 

PRF  PLOT 

PRFC 

ON 

NPRF 

7 

NIPRF 

6 

PRF  DATE 

PRFD 

158.8 

PRFD 

172.8 

PRFD 

200.8 

PRFD 

214.8 

PRFD 

228.8 

PRFD 

264.8 

PRFD 

284.8 

PRFD 

PRFD 

PRF  FREQ 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

100.0 

PRFF 

PRFF 

PRF  SEG 

IPRF 

15 

IPRF 

18 

IPRF 

22 

IPRF 

26 

IPRF 

28 

IPRF 

31 

IPRF 

IPRF 

IPRF 

TSR  PLOT 

TSRC 

ON 

NTSR 

1 

TSR  DATE 

TSRD 

144.5 

TSRD 

TSRD 

TSRD 

TSRD 

TSRD 

TSRD 

TSRD 

TSRD 

TSR  FREQ 

TSRF 

0.25 

TSRF 

TSRF 

TSRF 

TSRF 

TSRF 

TSRF 

TSRF 

TSRF 

VPL  PLOT 

VPLC 

OFF 

NVPL 

7 

VPL  DATE 

VPLD 

144.5 

VPLD 

158.5 

VPLD 

172.5 

VPLD 

200.5 

VPLD 

214.5 

VPLD 

228.5 

VPLD 

264.5 

VPLD 

284.5 

VPLD 

VPL  FREQ 

VPLF 

100.0 

VPLF 

100.0 

VPLF 

100.0 

VPLF 

100.0 

VPLF 

100.0 

VPLF 

100.0 

VPLF 

100.0 

VPLF 

100.0 

VPLF 

CPL  PLOT 

CPLC 

OFF 

NCPL 

1 

CPL  DATE 

CPLD 

212.208 

CPLD 

CPLD 

CPLD 

CPLD 

CPLD 

CPLD 

CPLD 

CPLD 

CPL  FREQ 

CPLF 

1.0 

CPLF 

CPLF 

CPLF 

CPLF 

CPLF 

CPLF 

CPLF 

CPLF 

RESTART 

RSOC 

OFF 

NRSO 

1 

RSIC 

OFF 

RSO  DATE 

RSOD 

172.5 

RSOD 

RSOD 

RSOD 

RSOD 

RSOD 

RSOD 

RSOD 

RSOD 

RSO  FREQ 

RSOF 

100.0 

RSOF 

RSOF 

RSOF 

RSOF 

RSOF 

RSOF 

RSOF 

RSOF 

CST  COMP 

CCOKPC 

ON 

LIMC 

OFF 

SDC 

ON 

FREQUK 

1 
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All 


CST  ACT 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

ACC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CST  ICON 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

CIC 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CST  PENT 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

CPRNC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CIN  CON 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

INACC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CTR  CON 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

TRACC 

OFF 

OFF 

OFF 

ON 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

COT  CON 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

DTACC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

CPR  CON 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

PRACC 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

OFF 

EX  COEF 

EXH20 

EXINOR 

EXORG 

BETA 

0.45 

0.01 

0.3 

0.45 

CO LI FORM 

COLQIO 

COLDK 

1.04 

1.4 

S  SOLIDS 

SSETL 

2.0 

ALGAE 

AGROW 

AMORT 

AEXCR 

ARESP 

ASETL 

ASATUR 

ALGDET 

1.5 

0.05 

0.02 

0.02 

0.14 

50.0 

0.80 

ALG  RATE 

AGT1 

AGT2 

AGT3 

AGT4 

AGK1 

AGK2 

AGK3 

AGK4 

10.0 

30.0 

35.0 

40.0 

0.1 

0.98 

0.98 

0.1 

DISS  ORG 

LABDK 

LRFDK 

REFDK 

0.12 

0.001 

0.001 

DETRITUS 

DETDK 

DSETL 

0.08 

0.35 

ORG  RATE 

OMT1 

OMT2 

OMK1 

OMK2 

4.0 

20.0 

0.1 

0.98 

SEDIMENT  SEDDK 
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S  DEMAND 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

SOD 

0.050 

0.050 

0.050 

0.10 

0.150 

0.20 

0.250 

0.250 

0.250 

0.250 

0.30 

0.350 

0.350 

0.350 

0.350 

0.350 

0.350 

0.30 

0.250 

0.250 

0.250 

0.250 

0.250 

0.20 

0.20 

0.20 

0.20 

0.150 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.150 

0.150 

0.10 

0.050 

WOOD  TEMP 

TWCOD 

1.0147 

WC  DEMAND 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

WCOD 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.70 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.40 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

0.30 

CBOD 

KBOD 

TBOD 

RBOD 

0.25 

1.0147 

1.85 

PHOSPHOR 

P04REL 

PARTP 

AHSP 

0.015 

1.2 

0.009 

AMMONIA 

NH3REL 

NH3DK 

PARTN 

AHSN 

0.08 

0.12 

1.0 

0.014 

NH3  RATE 

NH3T1 

NH3T2 

NH3K1 

NH3K2 

5.0 

20.0 

0.1 

0.98 

NITRATE 

N03DK 

0.12 

N03  RATE 

N03T1 

N03T2 

N03K1 

N03K2 

5.0 

20.0 

0.1 

0.98 

SED  CO 2 

C02REL 

0.1 

IRON 

FEREL 

FESETL 

0.5 

2.0 

STOICHMT 

02NH3 

020RG 

02RESP 

02ALG 

BIOP 

BION 

BIOC 

4.57 

1.4 

1.1 

1.4 

0.011 

0.08 

0.45 

02  LIMIT  02LIM 

0.00 


BTH 

FILE . 

. BTHFN 

bth.npt 

VPR 

FILE . 

. VPRFN 

vpr.npt 

LPR 

FILE . 

lpr.npt 

RSI 

FILE . 

rsi.npt 
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MET  FILE 


METFN 


set. npt 

QWD  FILE . 

not  used 

QIN  FILE . 

Br  1  qin_.br  1 .  npt 

Br  2  qinj»r2.npt 

TIN  FILE. . . 

Br  1  tin_.br  1.  npt 

Br  2  tin_br2 .npt 

CIN  FILE . 

Br  1  cin_brl.npt 

Br  2  cinjbr2.npt 

QOT  PILE . 

Br  I  qot_.br  1.  npt 

Br  2  qot_.br  2.  npt 

QTR  FILE . 

Tr  1  qtr_.tr  1 .  npt 

TTR  FILE . 

Tr  l  ttr_trl.npt 

CTR  FILE . 

Tr  1  ctr_trl.npt 

QOT  FILE . 

Br  1  not  used 

Br  2  not  usod 

TOT  FILE . 

Br  X  not  usod 

Br  2  not  used 

CDT  FILE . 

Br  X  not  used 

Br  2  not  usod 

PRE  FILE . 

Br  1  not  usod 

Br  2  not  usod 

TPR  FILE . 

Br  X  not  usod 

Br  2  not  usod 

CPR  FILE . 

Br  1  not  usod 

Br  2  not  usod 

EUH  FILE . 

Br  X  not  used 

Br  2  not  used 


QWDFH, 

QINFN. 

TINFNi 

CINFN. 

QOTFN. 

QTRFN 

TTRFN. 

CTRFN 

QDTFN 

TDTFN 

CDTFN 

PREFN 

,  TPRFN 

.  CPRFN 

.  EUHFN 
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nnt  ptt.*  . 

. TUHFH . 

Br  l  not  used 

Br  2  not  used 

CVH  FILE . . 

. . . CUHFN . . . 

Br  1  not  used 

Br  2  not  used 

BOH  FILE . 

Br  1  not  used 

Br  2  not  used 

TDH  FILE . 

Br  1  not  used 

Br  2  not  used 

mu  FTIiK.  . . . 

. CDHFH . . . 

Br  1  not  used 

Br  2  not  used 

SNP  FILE . . . 

. SNPFN . 

snp_run29 . opt 

TSR  FILE . 

tsr_run29 . opt 

PRF  FILE . 

prf_run29 . opt 

VPL  FILE . 

vpl.opt 

rpi.  fii* . 

. CPLFN . 

cpl.opt 
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Figure  B1 .  Sensitivity  analysis  results  from  increasing  SOD  parameter  50  per¬ 
cent  for  1981  (Sheet  1  of  5) 
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Figure  B2.  Sensitivity  analysis  results  from  decreasing  SOD  parameter  50  per¬ 
cent  for  1981  (Sheet  1  of  5) 
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Figure  Cl .  Scenario  1  results  from  increasing  pool  1 1  ft  for  1981  and  1983 
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